4 He adsorbed inside (10,10) single walled carbon nanotubes 



M.C.Gordillo 

Departamento de Sistemas Fisicos, Quimicos y Naturales. Facultad de Ciencias Experimentales. 
Universidad Pablo de Olavide. Carretera de Utrera, km 1. 41013 Sevilla. Spain. 



o 
o 

(N 
> 

o 

00 



J. Boronat and J. Casulleras 
Departament de Fisica i Enginyeria Nuclear, Universitat Politecnica de Catalunya, 
B4-B5 Campus Nord. 08034 Barcelona. Spain. 
(Dated: February 2, 2008) 

Diffusion Monte Carlo calculations on the adsorption of 4 He in open-ended single walled (10,10) 
nanotubes are presented. We have found a first order phase transition separating a low density 
liquid phase in which all 4 He atoms are adsorbed close to the tube wall and a high density arrange- 
ment characterized by two helium concentric layers. The energy correction due to the presence of 
neighboring tubes in a bundle has also been calculated, finding it negligible in the density range 
considered. 

PACS numbers: 68.90.+g,05.30.-d 
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One of the most interesting capabilities of carbon nan- 
otubes and their bundles is their purported use as reser- 
voirs, in particular for the lightest quantum fluids ( 4 He, 
H2 and Ne)ii2i2i£i5i£. The most obvious adsorption places 
are the outer part of the bundles, formed by the depres- 
sions known as grooves between every two tubes, and 
the rest of the exposed external surfaces. These ad- 
sorption places have already been experimentally doc- 
umented 3 ^^. One can think also of the interchannels 
located every three - or more if the arrangement is irreg- 
ular - tubes in the inner part of the bundle as accessible 
sites, but this latter possibility seems uncertain since any 
single vacancy in the carbon structure of the nanotubes 
could block those very thin interstices^. 

The other possible adsorption site is the inner part of 
the tubes themselves. This would imply to remove the 
ending caps of the nanotubes and thus allowing the gases 
to enter. To our knowledge, there is no experimental in- 
formation for this last possibility in the case of 4 He, the 
majority of theoretical studies being related to the limit 
of infinite dilution 1 ^ or to quasi one-dimensional narrow 
environment o 11 ! 13 ! 14 , even though there are some excep- 
tion a 15 ' 16 : 17 . In the present work, we present a study of 
4 He adsorbed in a (10,10) single walled carbon nanotube, 
one of the most common types. We have also estimated 
the effect that neighboring tubes have in the adsorption 
energies of 4 He inside a single tube. That can help to 
make a good estimation of the real adsorption capabili- 
ties of nanotube bundles. 

We have carried out a microscopic study of the sys- 
tem at zero temperature using the Diffusion Monte Carlo 
(DMC) method. DMC is a stochastic technique for solv- 
ing the many-body Schrodinger equation, providing for 
boson systems exact results within some statistical uncer- 
tainties 1 ^. The carbon atoms in the nanotubes have been 
explicitly considered, meaning that the carbon nanotube- 
hclium potential was obtained by summing up all the 
C-He pair interactions (taken from Ref. Ilfj), with the 
carbon atoms in fixed positions. This is not the most 



common setup, since most of the previous studies were 
carried out considering a smoothed helium-tube inter- 
action (an exception being Ref. Il5r ). The helium-helium 
potential is the HFD-B(HE) Aziz potential^. For the liq- 
uid phase simulations (see below), we have used N = 150 
4 He atoms, and to run calculations at different densities 
we have changed the length of the simulation cell between 
1537.5 and 34.44 A, these values being in all cases a mul- 
tiple of the elementary carbon cell of length 2.46 A, char- 
acteristic of an armchair nanotube. When a solid-liquid 
system was considered, we fixed the number of atoms 
belonging to the solid phase, and varied the number of 
atoms in the liquid one. Periodic boundary conditions 
in the transverse direction (z in cylindrical coordinates) 
have always been assumed. 

The trial wave function needed in the method for guid- 
ing the diffusion process was chosen to be of the form: 
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with N and N c the number of 4 He and C atoms, re- 
spectively, and ri the radial distance of particle i to the 
center of the nanotube. The first two terms of VP in Eq. 
(fTj) correspond to two-body correlation factors which ac- 
count for dynamical correlations induced by the He-He 
and He-C potentials. Both f He-Heir) and /hc-cM are 
of McMillan type, /(r) = exp[-0.5(6/r) 5 ]. The value 
for the parameters bue-He has been taken form Ref. [TR 
The value of &c-Hc varied from the phases with one he- 
lium layer (&Hc-Hc = 0) to the ones with two helium 
layers (6hc-Hc = 2-3 A, see below). In both cases, b was 
obtained by means of separate Variational Monte Carlo 
(VMC) calculations and found to be independent on the 
helium density. 

The variance of the calculation is reduced by introduc- 
ing also in the trial wave function (TT]) a one-body term 
3>(fi) which accounts for the mean interaction between a 
single He atom and the nanotube wall. At small densities, 



2 



we have taken for <E>(ri) the solution of the Schrodinger 
equation for a single 4 He atom in the cylindrically av- 
eraged potential of Ref. HH for a tube of the same di- 
ameter than a (10,10) one. This function has a single 
peak located roughly in the minimum of the considered 
helium-tube potential. When the density increases one 
can see the emergence of a second peak in the density 
profile obtained from the DMC simulations. This peak 
is located nearer to the center of the nanotube than the 
previous one. Thus, in the high density regime we veri- 
fied that the variance is further reduced by changing the 
low-density one-body term to another one which contains 
this closer-to-the-center peak in addition to the one near 
the carbon wall. In particular, we have used for <E>(ri) 
the square root of the radial density profile obtained in 
a DMC calculation at density p = 0.028 A" 3 . This trial 
wave function has two maxima and a reduced density 
(but not zero) between them. In principle, this profile 
with the double peak would correspond to a liquid phase 
since no localization on fixed sites of any of the helium 
atoms was imposed. 

When the density increases even more the layer closer 
to the wall can crystallize. In order to characterize this 
possible solid layer we have also studied the system by 
introducing Nosanow- Jastrow terms describing a solid or- 
der around the cylinder, 

* 8 (r u ...,r N ) = JJffflri-Ril) *(n, . . . , r N ) , (2) 

i=l 

with 'F the trial wave function of the liquid phase with 
a double peak (JTJ and g(r) a gaussian factor localizing a 
given set of particles around predefined lattice sites R^. 
The lattice chosen for the simulation corresponds to a 
triangular-like planar lattice rolled up to form a cylinder 
concentric to the carbon nanotube. A VMC optimization 
process was performed first by using a fixed number of 
atoms (150) for different lengths of the simulation cell, 
and varying also the number of atoms located in the first 
and second helium layers. We then compared the simu- 
lation results to those without the gaussian constrains at 
the same density. Those set of optimizations leaded us to 
conclude that the minimum energy configurations in the 
density range studied correspond to a lattice formed by 
slices of eight atoms that approximately share the same z 
coordinate (that of the tube axis) and are separated by a 
distance of 2.46 A in the z direction. This layer is located 
at a distance of 2.8 A from the carbon wall. Any other 
2D solid phases, including registered ones we found to 
produce configurations of higher energy. In the final sim- 
ulations, we have used 14 nanotube cells that, with this 
optimized structure, represent a solid layer composed by 
112 4 He atoms. The rest of the helium, up to 150 atoms, 
was considered to be in the second layer and was disposed 
without gaussian restrictions. The free parameter of the 
gaussian g(r) = cxp(— or 2 ) has been optimized (a=2.30 
A~ 2 ) showing negligible dependence on the density in the 
range studied. A check was also made to see the influ- 
ence of the corrugation imposed by the carbon structure 
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FIG. 1: Energy versus density for different elections of the 
trial function: One layer liquid (full circles); (squares); two 
layer liquid (circles); a close to the wall solid (full squares). 



in the energy of the solid layer. To do so, we located the 
solid lattice in different positions relative to the underly- 
ing structure. The energies obtained were similar within 
the error bars, what indicates that the solid layer is a 
two-dimensional incommensurate solid. 

The results for the total set of simulations of 4 He in- 
side a single (10,10) tube are displayed in Fig. 1. The 
densities were calculated by straightforwardly dividing 
the number of atoms by the nanotube volume, consid- 
ering a radius of 6.8 A in all cases. Full circles are the 
results for a one-layer He phase, open circles correspond 
to the two-layer liquid arrangement, and the full squares 
indicate the energies obtained for the phase with a solid 
layer closer to the carbon wall. The main conclusion one 
can draw from this figure is that there is a minimum of 
energy of the one-shell phase; corresponding to zero pres- 
sure and that can be found from a third order polynomial 
fit to be at 0.0086 A -3 with a corresponding binding en- 
ergy of 187.18 ± 0.02 K. The same binding energy for the 
infinite dilution limit is 185.01 ± 0.02 K. That is bigger 
than the experimental value for graphite, 143 ± 2 K 25 , 
and smaller than the corresponding values for a (5,5) 
tube (429.984 and 429.966 K at zero pressure and for in- 
finite dilution, respectively). This simply means that an 
increasing in the tube radius implies a lowering of the 
binding energy to approximate the one in the flat surface 
of graphite. What is more important, the stable phase 
at zero pressure and K is a liquid. In this respect, the 
situation is different than in the case of graphite 2 ^ and 
for the cylindrical system considered in Ref. Ilq . 
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FIG. 2: Maxwell construction to determine the stability range 
of the phases involved in the modeling of the system. Full line, 
one layer liquid and solid+liquid phases ; Circles, two layer 
liquid phase. 
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FIG. 3: Radial density profiles at zero pressure (dashed line), 
at the lower and upper limits of the coexistence regions (dot- 
ted and long dashed line, respectively). The solid line repre- 
sents the case for the solid+liquid phase at the highest density 
simulated. 



The stability of the rest of the phases can be deduced 
with the help of Fig. 2. This is a Maxwell construction 
for the density range in which we could have a two-layer 
liquid or a two-layer solid+liquid. Full lines are third- 
grade polynomial fits to the corresponding one layer liq- 
uid phase, and the two-layer solid+liquid arrangement, 
while circles are the results for the two-layer liquid. The x 
axis corresponds to the inverse of the density as defined in 
the previous figure. The dashed line is a double tangent 
construction between the one-layer liquid and the phase 
including a solid. In a conservative frame, this would 
indicate a transition between a low density phase of 2 
10~ 2 A" 3 and a one including a solid, of 2.7 10~ 2 A~ 3 . 
The transition pressure would be of 430 atm. However, 
another alternative explanation can be deduced from the 
data, and it is a double transition between a single layer 
liquid and a double layer liquid, followed by another one 
between this last dense liquid and the solid and liquid 
phase. The corresponding coexistence densities are the 
same for the two phases already mentioned, but we will 
have small a window with a two-layer liquid at around 2.5 
10~ 2 A~ 3 . The one shell liquid- two shell liquid transition 
would be then at 420 atm. Unfortunately, the available 
data (and shown in Table I) will not allow us to discrimi- 
nate between both scenarios. In any case, it is clear that 
an increasing in the density will eventually create a solid 
close to the wall. 

To complete the characterization of the phases, we dis- 
played the radial profiles for different helium densities 
in Fig. 3. The dashed line represents the case of zero 
pressure (0.0086 A -3 ), while the dotted line displays the 
situation in the low density limit of the transition region. 
One can see that in both cases the system is formed by a 
single cylindrical helium shell, located at around 4 A from 



the center of the tube. The main difference is the increas- 
ing in the helium density that can be inferred from the 
growth in the corresponding peak. The long-dashed line 
indicates the situation in the upper density limit of the 
coexistence zone for the case of the solid+liquid phase. 
We observe a second peak located closer to the center 
of the tube, and another one corresponding to the two- 
dimensional solid, roughly in the same position than in 
the single layer phase. Finally, the full line represents 
the solid+liquid phase at the highest density simulated. 
There, the first peak is virtually identical to the one at 
the transition, as befits to a solid of unchanging density, 
going the additional helium to fill the second layer. 

Until now, we considered only the case of an isolated 
carbon nanotube, but this is not the most common situ- 
ation, since those cylinders tend to group to form what is 
termed a bundle. There, each tube is ideally surrounded 
by six others, forming a structure whose section is sim- 
ilar to a triangular 2D lattice with a minimum distance 
between centers of cylinders of 17 A3 2 -. Taking that into 
account, we estimated the influence that the surrounding 
tubes have in the energy per helium atom within a mean 
field scheme by means of the following expression: 

Ecorrection = 6/ / d(x' , y')dx' dy' (3) 

J x' J y' 

x / / R(x,y,z)V(x,y,z,x',y')dxdydz 

JO J x Jy 

where d(x',y') represents the normalized probability of 
finding an 4 He atom at coordinates x' and y' for any z 
position in the first tube. R{x, y, z) is the corresponding 
radial density function as the ones depicted in Fig. 3, 
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TABLE I: Energies per 4 He atom in K at several densities. 



Density (A~ 3 ) 


(10,10) tube 


E correction 


Low density phase 


6.7 10" 4 


-185.02(4) 


-0.003 


1.7 1CT 3 


-185.27(3) 


-0.006 


3.4 1(T 2 


-185.81(3) 


-0.013 


4.2 1CT 3 


-186.14(4) 


-0.016 


5.6 1CT 3 


-186.58(3) 


-0.022 


7.0 10~ 3 


-186.93(4) 


-0.027 


8.4 10~ 3 


-187.14(3) 


-0.033 


1.05 10~ 2 


-186.89(3) 


-0.041 


1.2 10~ 2 


-186.02(4) 


-0.047 


1.4 10" 2 


-183.76(4) 


-0.055 


1.6 10~ 2 


-180.40(4) 


-0.062 


1.7 10" 2 


-176.66(4) 


-0.067 


1.8 10~ 2 


-169.97(4) 


-0.073 


1.9 10" 2 


-165.59(6) 


-0.077 


2.0 10~ 2 


-159.45(7) 


-0.081 


Two layer liquid density phase 


1.8 10 -2 


-155.60(9) 


-0.069 


1.9 10~ 2 


-152.56(9) 


-0.072 


2.0 10~ 2 


-149.15(7) 


-0.075 


2.1 10~ 2 


-146.22(9) 


-0.079 


2.2 10~ 2 


-141.58(6) 


-0.083 


2.3 10 -2 


-136.66(6) 


-0.087 


2.5 10~ 2 


-130.22(8) 


-0.092 


2.6 10~ 2 


-122.49(8) 


-0.099 


2.8 10~ 2 


-112.2(1) 


-0.103 


Two layer liquid+solid density phase 


2.2 10-' J 


-134.58(5) 


-0.092 


2.4 10" 2 


-129.63(5) 


-0.095 


2.6 10~ 2 


-124.50(5) 


-0.099 


2.7 10~ 2 


-118.30(5) 


-0.103 


2.9 10~ 2 


-109.64(6) 


-0.107 


3.0 10~ 2 


-98.41(9) 


-0.111 



and depends on the location of helium atoms in the sec- 
ond tube {x, y, z). Finally, V(x, y, z, x' , y') is the helium- 



helium Aziz potential used in our simulations? . The 
result of the integral is multiplied by six to take into ac- 
count the number of tubes surrounding a particular one 
at the same distance. A check of the accuracy of this 
kind of mean field estimation has already been done for 
a closely related system^ 4 an d found a very good approx- 
imation to the results of a full Monte Carlo calculation. 

The results of this energy correction are given in Table 
I. One thing is immediately apparent, and it is that the 
mean field term increases with density, but in all cases is 
negligible and of the order of the error bars of the DMC 
calculation. That means that the binding energy of 4 He 
in a bundle is essentially the same that the one for a 
single tube. 

Summarizing, we have performed a full DMC calcu- 
lation on the subject of 4 He adsorbed in carbon nan- 
otube bundles. We have found a low-density liquid phase 
characterized by a single helium layer close to the nan- 
otube wall. When the density increases there is a phase 
transition to a two-layer arrangement that will eventu- 
ally solidify with further increase of the number of ad- 
sorbed helium atoms. We also saw that the influence 
of the neighbouring tubes in a bundle in the properties 
of a given one is negligible. Experimental signatures of 
the phases predicted in this work could be obtained by 
means of low temperature specific heat measurements 7 
and neutron diffraction 6 . 
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